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ABSTRACT 



O 
O 

(N , 

j>^! We present SMARTS consortium optical/IR light curves of SN 2006aj asso- 

ciated with GRB 060218. We find that this event is broadly similar to two pre- 
viously observed events SN 1998bw/GRB 980425 and SN 20031w/GRB 031203. 
^ ! In particular, all of these events are greatly under-luminous in gamma-rays com- 

pared to typical long-duration GRBs. We find that the observation by Swift of 
even one such event implies a large enough true event rate to create difficul- 



CN ! ties in interpreting these events as typical GRBs observed off-axis. Thus these 

QQ i events appear to be intrinsically different from and much more common than 

^ i high-luminosity GRBs, which have been observed in large numbers out to a red- 

\^ I shift of at least 6.3. The existence of a range of intrinsic energies of GRBs may 



present challenges to using GRBs as standard candles. 

Subject headings: gamma rays: bursts — supernovae: general — supernovae: 
individual (SN 2006aj) 



1. Introduction 

While some long-duration gamma-ray bursts (GRBs) are clearly associated with su- 
pernovae (SNe), a deeper understanding of the GRB/SN connection remains elusive. The 
GRB/SN link was first confirmed observationally with the detection of the low-reshift GRB 
980425/SN 1998bw {z = 0.0085) (Galama et al. 1998). GRB 980425, however, was not a 
typical GRB; it was under- luminous in gamma-rays and had no detected optical afterglow 
(OAG). SNe were later associated with typical GRBs at cosmological redshifts (e.g. Bloom et 
al. 1999; Delia Valle et al. 2003). However, another low gamma-ray luminosity event similar 
to GRB 980425/SN 1998bw was not observed until GRB 031203. This burst was also orders 
of magnitude under-energetic and, despite its low redshift {z = 0.1055) (Prochaska et al. 
2004), was followed by only a dim OAG (Malesani et al. 2004). Follow-up observations of 
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this burst detected a SN-like brightening (Cobb et al. 2004; Gal- Yam et al. 2004; Thomsen 
et al. 2004) and SN 20031w was confirmed spectroscopically by Tagliaferri et al. (2004). The 
spectra of SN 20031w were reminiscent of those of SN 1998bw (Malesani et al. 2004). The 
The two SN also had similiar peak magnitudes, although SN 20031w was somewhat brighter 
and evolved more slowly. Their light curve shapes were also qualitatively different, with SN 
1998bw climbing smoothly to peak while SN 20031w experienced a broad plateau. 

The low gamma-ray luminosity of GRB 980425 and 031203 suggested that they might 
represent a new GRB category. Alternatively, they could be normal GRBs that appear 
under-luminous because their jetted emission is observed off-axis (e.g. Yamazaki et al. 2003; 
Ramirez- Ruiz et al. 2005). A comparison of the event rate of low- to typical- luminosity 
GRBs was warranted but only a small and inhomogeneous sample of well-locahzed GRBs 
existed before Swift, as pice-Swift GRBs were detected using multiple instruments, each with 
unique sensitivity and sky coverage. Swift provides a large and homogeneous GRB sample 
that is well-suited for rate calculations. 

On 2006 February 18 at 03:34:30 UT Swift detected a new low-luminosity event: GRB 
060218 (Cusumano et al. 2006). This was an unusual GRB with weak gamma-ray emission 
lasting over 2000 seconds (Barthelmy et al. 2006). This burst was followed by an unusual 

OAG that brightened for 10 hours before decaying like a typical OAG (e.g. Marshall et al. 
2006). This was the first Swift GRB to be associated with a SN: SN 2006aj. The SN was 
initially noted in spectral observations (Masetti et al. 2006) and then detected as an optical 
re-brightening (e.g. D'Avanzo et al. 2006; Ovaldscn et al. 2006, etc.). At z = 0.033, GRB 
060218 is now the second closest GRB with a measured redshift (Mirabal & Halpern 2006). 
This is the third example of a GRB-related SN in which the gamma-rays are highly under- 
luminous and the SN fight curve is clearly distinct from the GRB's OAG. Hereafter, we will 
refer to these Long-duration, Low- Luminosity events as L3-GRBS. 

SMARTS observations of SN 2006aj began on 2006 Feb 22 at 00:35 UT (Cobb k Bailyn 
2006). We present optical/IR data obtained between 5 and 30 days following GRB 060218. 
Our homogeneous data demonstrate that the light curve of SN 2006aj is quafitatively similar 
to that of the pre-Swift L3-GRB SNe 1998bw/20031w. We argue in §4 that the detection of 
this single event in the Swift era already places strong constraints on the nature of L3-GRBS. 
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2. Observations and Data Reduction 

Our data was obtained using the ANDICAM instrument mounted on the 1.3m telescope 
at Cerro Tololo Inter-American Observatory.^ This telescope is operated as part of the 
Small and Moderate Aperture Research Telescope System (SMARTS) consortium.^ Nightly 
imaging was obtained over 26 days with occasional interruptions for weather and equipment 
problems. The GRB/SN was only observable for a limited period of time immediately after 
twihght (< Ihr). Consequently, all observations were obtained at high airmass (sec (z) > 2). 

Each nightly data set consisted of 6 individual 360-second I-band observations obtained 
simultaneously with 30 dithered 60-second J-band images. The data were reduced in the 
same way as in Cobb et al. (2004). A few additional steps were added, including cosmic ray 
removal in the I-band images using the L.A. Cosmic program^ (van Dokkum 2001) and I- 

band fringe correction using an iterative masking technique. Some images were not included 
in the final frames because of excessive background due to twilight or because of telescope 
drift. Typically, the final frames were equivalent to 30 minutes of 1/ J-band exposure time. 

The relative magnitude of the SN -|- host galaxy was determined by comparison with 11 
(3) on-chip, non-variable objects in I (J) using seeing matched aperture photometry. Differ- 
ential magnitudes were converted to apparent magnitudes by comparison, on photometric 
nights, with Landolt standard stars in the fields of RU149 and PG1047 (Landolt 1992) for 
the I-band images, and with 3 on-chip 2MASS stars (Skrutskie et al. 2006) for the J-band 
images. The difference in airmass value between the Landolt standard frames and science 
frames was corrected for using an extinction coefficient of 0.066 magnitudes per airmass. 

The light curves are shown in Figure 1 and the photometric data are summarized in 
Table 1. The error bars represent the photometric measurement error, which accurately 
reflects nightly variations in image quality but does not account for systematic measurement 
errors. In addition to the relative night-to-night uncertainty, there is a systematic error of 
0.05 magnitudes in I and J resulting from uncertainties in the photometric calibration. 



^http: / /www.astronomy.ohio-state.edu/ANDICAM 

^http : / /www. astro .yale.edu /smarts 

^http: / / www.astro.yale.edu/dokkum/lacosmic/ 
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3. Results 

Figure 1 shows that GRB 060218's optical and IR counterpart brightened for the first 
two weeks and then proceeded to gradually decay. This behavior is not consistent with that 
of a standard GRB OAG (e.g. Tagliaferri et al. 2005) but is reminiscent of the low redshift 
events SNe 1998bw/20031w. Identification of this optical emission as a SN is possible from 
our data alone due to our dense observations and the object's particular transient behavior; 
spectral evidence obtained by other groups clearly supports our claim (Modjaz et al. 2006; 
SoUerman et al. 2006; Mirabal et al. 2006). The well-sampled nature of our observations 
allows us to determine an unambiguous time of peak brightness in I and J. This parameter 
will be important for determining the amount of mass ejected in the supernova explosion, 
thought this modeling is beyond the scope of this paper. The position of peak brightness 
was determined by fitting second order cubic splines to the data points, with errors derived 
from the formal chi-squared error on the fits in combination with the error on the measured 
magnitudes. The combination of the host galaxy and SN reaches a peak apparent magnitude 
in I of 16.91 ± 0.05 mag after 13.1^11 days and in J of 16.65 ± 0.06 mag after 17.6l|^ days. 
The rest-frame time to peak is, therefore, approximately 12.7 days in I and 17.0 days in J. 

The Galactic extinction correction along the line of sight to the host galaxy is taken to 
be Aj — 0.23 mag and Aj — 0.11 mag, assuming the Galactic extinction curves of CardeUi, 
Clayton, & Mathis (1989) and a measured reddening value of E(B-V)=0.127 mag (Guenther 
et al. 2006). The pre-burst SDSS model magnitude of the host galaxy is i = 19.805 ± 0.041 
mag, not corrected for Galactic extinction (Adelman-McCarthy et al. 2006). Using the 
transformation equations derived by Lupton (2005), this corresponds to / = 19.368 ± 0.047 
mag. The peak absolute magnitude of the supernova is, therefore, Mj = —19.02 ±0.09 mag. 
No k-correction has been applied, but this should result in minimal error due to the low 
redshift of the burst. A correction of -0.04 magnitudes was applied to account for spectral 
stretching. The exact pre-burst J magnitude of the host is unknown as the host galaxy is 
too dim to appear in the 2MASS catalog. Our observations indicate the host galaxy must 
have a J magnitude > 18. Assuming a range of host magnitudes from 18 mag to 20 mag, 
the peak absolute magnitude of the supernova in J is approximately Mj — —19.1 ±0.2 mag. 

Note that 2006aj clearly peaks later in J than in I. This later peak at redder wavelengths 
follows the trend seen in SN 1998bw, which, in the rest-frame, peaked 1.6 days earlier in 
V than in I. The rest-frame V-band peak of SN 2006aj occurred at approximately 9.7 days 
(Modjaz et al. 2006), which is 2.9 days prior to the I-band peak. Likewise, SN 20031w peaked 
in V at ~ 18 days (Malesani et al. 2004) and in I at ~ 23 days (Cobb et al. 2004; Malesani et 
al. 2004). The combined light of the galaxy and the SN reddens from I- J = 0.0 mag during 
the first week to I-J = 0.6 mag for the last few observations. This is a stronger evolution 
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in I-J color than experienced by either SN 1998bw or SN 20031w, whose I-J colors in the 
first month only change by about 0.3 magnitudes (Gal- Yam et al. 2004). This comparison 
is complicated, however, by the unknown intrinsic I-J color of the host galaxy of SN 2006aj. 



Our data, together with those of Modjaz et al. (2006), SoUerman et al. (2006) and 
Mirabal et al. (2006), show that GRB 060218 is the third detected GRB to be followed by 
a dominant SN. With such a small sample being used to extrapolate the characteristics for 
an entire category, it is important to collect homogeneous and detailed observations over a 
wide range of wavelengths. Our data provides dense IR coverage, which extends the total SN 
2006aj dataset out to redder bands than reported thus far. It is instructive to compare all 
three cases in which Ls-GRBs have been detected, each of which was followed by a Type Ic 
SN: 1998bw, 20031w and 2006aj (see Figure 2). We note that the hmit on observing L3-GRB 
events are more stringent than those on observing the associated SNe, so it is unlikely that 
GRBs for which no optical counterpart are observed are of this character. The properties of 
these three events are shown in Table 2. All three SNe are very similar in peak brightness, 
though 20031w may be half a magnitude brighter than the others. The biggest difference 
between the bursts is their rise times, with 2006aj peaking the fastest and 20031w taking 
the longest time to peak. The rest-frame photon energy at which the GRB spectrum peaks 
(Ep^i) appears to increase with increasing SN rise time. 

As discussed in the introduction, it would be instructive to compare the frequency of 
these L3 bursts with that of high-luminosity bursts using the homogeneous Swift dataset. 
Such a comparison is now possible since GRB 060218 is the first ^m/t- detected burst that 
falls in the category of low-luminosity GRBs associated with Type Ic SNe. The low observed 
fiuxes and redshifts of L3-GRBS suggest that the underlying rate of L3-GRB events may be 
quite high (see also Pian et al. 2006; Soderberg et al. 2004), since the volume in which these 
sources can be observed is much smaller than that of typical GRBs. The ratio of the event 
rate of L3-GRBS to ordinary long-duration GRBs is expected to be 



where Ri^t denotes the true rate per comoving volume of the two kinds of events, Robs is the 
observed event rate seen by a given experiment, and is the comoving radial distance out 
to which the events could be observed by that experiment. In the case of Swift, the limiting 
volume in which L3-GRBS are observable is so small that the observation of even one of 
these events implies an extremely high ratio of true rates. 



4. 



Discussion 
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Specifically, GRB 060218 would not have triggered the BAT event monitor if it had been 
~ 2 times fainter, which corresponds to a maximum redshift of ^ = 0.046. By contrast, the 
31 high-luminosity, long-duration Swift GRBs that have measured redshifts have an average 
redshift of z = 2.6. Many of these bursts would have been detected even if they had occurred 
ai z > 6. We will take a conservative approach, however, and assume that most of these 
bursts would not have been detected had they been at extreme redshifts and adopt z = 2.6 
as the redshift below which the high-luminosity GRB samples are complete. Assuming the 
concordance cosmology of Qa — 0.73, Qm — 0.27 and a Hubble constant of 71 km s~^ Mpc~^, 
the ratio of the comoving radial distance cubed for these two events is 3.1 x 10^. Assuming 
that the ratio of rates is equal to the number of long-duration GRB events observed to 
date (1 L3-GRB per about 100 GRBs), we find a true event ratio of 3 x 10^. In the near 
future, thanks to Swift, the cosmological GRB sample may be complete out to at least a 
redshift of 2; ~ 5. For this volume, the true event ratio increases by a factor of ~ 2. If the 
current highest GRB redshift measurement, z = 6.3, is representative of the redshift out to 
which we have a complete sample, then the true event ratio is ~ 10^. We also note that we 
are assuming a constant GRB rate per unit time. If, instead, the rate scales with the star 
formation rate (e.g. Natarajan et al. 2005), then the relevant ratio would be higher still. 

These event ratios provide a constraint on the origins of L3-GRBS. One explanation of 
these events is that they are standard GRBs observed off-axis (e.g. Yamazaki et al. 2003; 
Ramirez- Ruiz et al. 2005). This is an attractive option as it accounts for all long-duration 
GRBs using a "unified model" , with observed differences attributed only to viewing angle. 
However, this scenario implies a maximum true rate ratio, which would be generated if the 
L3-GRBS could be observed from any angle. This upper limit is 27r/^|, where ^| is the jet 
opening solid angle of a typical GRB (in steradians). Jet breaks observed in X-ray/optical 
AG light curves constrain the jet opening angle to ~ 10° (e.g. Frail et al. 2001), so we infer 
a maximum true rate ratio of 65 in this model. This ratio is a factor of five lower than the 
ratio of 3 X 10^, which we inferred for z = 2.6, and lower by ~ 20 than for z = 6.3. 

Since Swift has observed only one L3-GRB, the underlying event rate remains uncertain; 
GRB 060128/SN 2006aj could have been a serendipitous event. From Poisson statistics, the 
90% lower confidence limit of the true event rate is an order of magnitude lower than assumed 
above, which implies that the off-axis scenario is still acceptable, provided one assumes a 
limit of z — 2.6 and a constant GRB rate. If Swift detected a second L3-GRB, the off-axis 
scenario would become significantly less probable. Such a Poissonian analysis addresses the 
question "given an event rate, what is the probability of seeing an event", whereas in this 
case one might more appropriately ask the Bayesian question "having seen an event, what 
is the probability of a given event rate". Assuming a uniform prior on the distribution of 
event rates (an assumption for which there is no real basis), we find that the off- axis scenario 
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is implausible at the 98% level. Of course, the undetermined luminosity function of GRBs 
could complicate this calculation, as could cosmic evolution of the GRB source population. 

At face value, however, the above calculation of the event rate suggests that there is 
a category of GRB events that is intrinsically different from that of typical GRBs. Several 
suggestions have been made for how these differences can be accounted for including the 

possibility that the gamma-rays arc produced in supernova shock breakout (Matzner & 
McKec 1999; Tan et al. 2001) or "failed collapsars" in which highly relativistic jets fail to 
develop due to baryon loading (Wooslcy & MacFadycn 1999). The orientation-corrected 
energies of GRBs have been claimed to be constant at ~ 10^^ ergs (Frail ct al. 2001). 
However, if intrinsically low-energy GRBs exist as a separate population, efforts to use GRBs 
as standard candles (e.g. Lazzati et al. 2006; Ghirlanda et al. 2004) may be compromised. 

We thank SMARTS observers D. Gonzalez and J. Espinoza for their dedication and S. 
Tourtellotte for assistance with optical data reduction. We greatly appreciate discussions 
with A. Cantrell and J. Emerson. This work is supported by NSF Graduate Fellowship 
DGE0202738 to BEG and NSF/AST grant 0407063 and Swift grant NNG05GM63G to GDB. 
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Table 1. Photometry of SN 2006aj in the Host Galaxy of GRB 060218 



Days after GRB^ 


I magnitude^ J magnitude^ 


4.88 


17.51 ±0.01 17.26 ±0.03 


5.87 


17.34 ±0.01 17.20 ±0.03 


6.87 


17.22 ±0.01 17.08 ±0.03 



^Days after burst trigger at 2006 Feb 18, 
03:34:30 UT. 

'^These values have not been corrected for 
Galactic extinction. There is an additional un- 
certainty of 0.05 mag in the transformation of 
relative to apparent magnitudes. 

Note. — The complete version of this table 
is in the electronic edition of the Journal. The 
printed edition contains only a sample. 
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Table 2. GRB/SN Properties 





GRB 980425 


GRB 031203 


GRB 060218 




SN 1998bw 


SN 20031w 


SN 2006aj 


redshift 


0.0085 


0.1055 


0.033 


flucnce (10~^ erg cm~^) 


2.8 ±0.5'' 


2.0 ±0.4'^ 


6.8 ±0.4^= 


total duration (s) 


~ 40 


~ 40 


> 2000 


Ep,i (keV) 


55 ± 15^ 


158 ±51^ 


< 10^ 


Eiso (1 X 10™ erg)^ 


0.010 ± 0.002'^ 


1.0 ± 0.4*^ 


0.65 ±0.15^ 


I-band Tp^afc (days)g 


17.7 iO.S'' 


18 - 28' 


12.71?:° 


peak Mj (mag) 


-19.27 ± 0.05*^ 


-19.0 to -19.7' 


-19.02 ±0.09 


1-J SN color. ~T,,(„A.. / 


0.5^ 


- 0. f 


- 0.0'" 



^(40 - 700 keV) Plan et al. (2000) 
^(20 - 200 keV) Sazonov et al. (2004) 
^(15 - 150 keV) Sakamoto et al. (2006) 
'lAmati (2006) 
^Amati et al. (2006) 

^rest-frame 1 - 10,000 keV, assuming Hq = 70 km s~^ Mpc~^, fl\ = 0.7, 
Qm = 0.3 

^in the rest-frame 

''Galama et al. (1998) 

'exact value depends strongly on extinction assumptions, Cobb et al. (2004); 
Gal- Yam et al. (2004); Malesani et al. (2004); Thomsen et al. (2004) 

'at ~ 22 days post-burst, Patat et al. (2001) 

'"assuming a J-band host magnitude of 19 
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Fig. 1. — Top panel: The I-band (circles) and J-band (squares) aperture photometry hght 
curve of the host+SN of GRB 060218. Values have been corrected for a Galactic foreground 
extinction of Ay =0.39 mag. For clarity, the I-band points have been shifted by +0.1 magni- 
tudes. Error bars are photometric measurement errors and do not include possible systematic 
effects. The curves are fit with second order cubic splines. Bottom panel: I-J color evolution, 
the combined light is observed to redden with time. 
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Fig. 2. — Absolute magnitude light curves of SN 1998bw (circles), SN 20031w (sqaures) and 
SN 2006aj (triangles) in I (top panel) and J (bottom panel). The SN 20031w data has been 
binned in intervals of 5 days. The apparent J-band host galaxy magnitude of SN 2006aj 
is assumed to be 19 mag. SN 20031w may be shifted by ~ —0.5 magnitudes if a stronger 
line-of-sight extinction is assumed. 



